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ABSTRACT 


The  prerequisite  to  a  dynamic  analysis  of  an  oc*an 
platform  is  an  accurate  baseline  measurement  of  the  fundamental 
flexular  and  torsional  natural  frequencies.  As  an  element 
of  this  baseline,  the  author  uses  the  Runga-Kutta  numerical 
integration  method  to  determine  the  influence  of  a  non¬ 
linear  structural  element  on  the  measurement  of  the 
fundamental  flexural  natural  frequency  of  the  system.  The 
non-linear  phenomena  which  was  investigated  was  the  collision 
between  the  vertical  conductor  pipes  and  their  guide  collars 
on  the  platform.  The  analysis  proved  that  for  large  excita¬ 
tion  loads  the  response  is  non-linear,  but  that  there  is  an 
excitation  level  below  which  the  non-linearity  is  negligible.. 
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CHAPTER  1 

Problem  Definition  and  Model  Description 

The  size  of  an  ocean  platform  and  the  variation  of 
environmental  excitations  on  it  make  the  determination  of  the 
integrity  of  the  platform  a  vital  concern  to  the  ocean 
engineer.  A  method  of  integrity  certification  which  is  under 
development  utilizes  on-location  measurement  of  the  platform 
natural  frequencies  to  determine  if  a  structural  failure  has 
occurred.  This  technique  requires  an  initial , accurate 
baseline  determination  of  the  platform  natural  frequencies. 
These  baseline  measurements  will  be  compared  to  values 
obtained  at  later  dates  to  establish  if  any  changes  have 
occurred.  However,  before  a  shift  in  natural  frequency  can 
be  related  to  damage,  all  other  sources  of  change,  such  as 
live  platform  load  variations  and  marine  growth,  must  be 
evaluated.  In  this  paper,  the  author  addresses  the  impact 
of  a  non-linear  structural  dynamic  element  on  the  measurement 
of  the  fundamental  natural  frequency. 

For  example  purposes,  the  platform  to  be  analyzed  is  a 
deep  ocean  drilling  platform  with  a  centrally  located, 
unattached,  conductor  which  extends  from  the  platform  deck 
to  the  ocean  floor.  The  purpose  of  the  analysis  is  to 
determine  the  effect  of  the  non-linear  interaction  between 
the  platform  and  the  conductor  on  the  measured  natural 
frequency  of  Che  entire  system  in  the  fundamental  flexural 
mode.  In  order  to  limit  the  problem  to  a  manageable  size. 
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several  constraints  must  be  placed  on  the  system.  Principal; 
among  these  is  the  use  of  a  periodic  excitation  to  drive 
the  system.  This  constraint  simplifies  the  identification 
of  the  platform  response  frequencies.  Secondly,  the  initial 
separation  between  the  platform  and  the  conductor  is  set  to 
zero.  In  the  author's  opinion,  this  position  presents  a  worst 
case  condition,  and  it  is  also  the  most  likely  from  engineering 
considerations  of  the  conductor  installation.  The  excitation 
loads  were  varied  in  order  to  demonstrate  a  breakpoint 
between  the  linear  and  non-linear  behavior  of  the  structure. 

The  remaining  system  parameters  were  set  equal  to  specific 
values  based  on  engineering  judgement  and  the  best  information 
available.  The  greatest  uncertainty  in  parameter  estimation 
occurred  in  the  determination  of  the  conductor  parameters. 

A  more  definitive  investigation  of  these  quantities  should 
be  performed  before  further  system  analysis  is  conducted. 

The  platform  can  be  modeled  as  a  mass,  spring  and  damper 
system  which  surrounds  another  mass,  spring  and  damper  system 
modeling  the  conductor.  The  platform  mass  is  free  to  move 
along  its  horizontal  axis,  and  if  sufficient  excitation  is 
provided,  collide  with  the  conductor  mass.  This  collision 
causes  a  reduction  in  the  platform  movement  while  increasing 
the  motion  of  the  conductor.  The  figure  on  the  next  page 
illustrates  the  model. 


■1  pi at fora  mass 

kl  platform  spring  constant 

bl  platform  viscous  damping  coefficient 

a 2  conductor  mass 

k2  conductor  spring  constant 

b2  conductor  viscous  damping  coeffidient 

xl  platform  displacement  from  equilibrium  position 
x2  conductor  displacement  from  equilibrium  position 


FIGURE  1.1 


SYSTEM  SCHEMATIC  DIAGRAM 
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A  mathematical  solution  for  the  response  of  tht  model  is  derived 
in  the  following  manner: 

1.  Equation  of  motion  of  the  platfcrm  moving 
separately : 


+  kxxx  =  F  sin  cot, 

or,  xi  =  (1/iUj)  (F  sin  cot  -  b1xl  -  kxxx)  (1) 

where  x  is  relative  to  the  unforced  equilibrium  position  of  the  platform. 
2.  Equation  of  motion  of  the  conductor  moving 
separately : 

m2x2  +  b2x2  +  k2x2  *  0 

or,  x2  *  (l/m2)  (-b2x2  -  k2x2)  (2) 

where  x  is  relative  to  the  unforced  equilibrium  position  of  the  conductor. 
.  3.  Equation  of  motion  for  the  combined  platform 

and  conductor  masses: 

(mi  +  m2)xi  ♦  biXx  +  b2xx  ♦  kxxx  +■  k2x2  *  F  sin  cot 
or,  xi  =  (l/mi+nu)  (F  sin  cot  -  (bi+b2)xi  -  kixi  -  k2x2)  (3) 

In  equations! , 2and  3,  masses,  damping  constants  and 
spring  constants  are  denoted  by  m,  b  and  k  respectively. 
Subscript  1  designates  the  platform  and  subscript  2  denotes 
the  conductor.  Equations  1  and  2  completely  describe  the 

■anjjpeiue  oft  tin  ayseemi  when  the-  pi and-c cnduc-t-ox -ar e - 

in  motion  separately.  Equation  3  describes  tho  response 
after  a  collision  has  occurred.  However,  several  questions 
about  the  system  have  yet  to  be  answered.  Specifically, 
what  happens  at  the  time  of  'ollislon,  how  is  the  time  of 
collision  determined,  and  how  is  the  time  of  separation 
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determined?  The  answers  to  these  questions  provide  the 
remaining  equations  which  will  define  the  motion  of  the 
platform/conductor  system  at  all  times. 

To  determine  when  a  collision  occurs,  the  gaps  between 
the  platform  and  the  conductor  must  be  continuously  determined 
When  the  value  of  the  gap  distance  on  either  side  of  the 
conductor  goes  to  zero,  a  collision  has  occurred.  The 
following  equation  defines  this  feature  mathematically. 

4.  Total  gap  distance: 

gapj  ♦  gap2  *  do  -  di  (4) 

where  do  is  the  platform  aperture  diameter  and  di  is  the 
conductor  outside  diameter. 

5.  Gap  distance  at  any  time,  t: 

gap^t)  •  gapl(t-l)  +  (x2(t)  -  x2(t-l)^ -(Xj(t) 

-  XxCt-1))  (5a) 

gap2(t)  =*  do  -  di  -  gap! (t)  (5b) 

where  (t-1)  defines  the  time  increment  just  prior  to  the 

time,t,  being  calculated. 

At  the  time  of  collision,  energy  will  be  lost.  However, 
momentum  is  conserved.  Therefore  the  velocity  of  the  system 
can  be  determined  as  shown  below. 
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6.  Conservation  of  momentum: 
mix  i  *  (m !  +  m2  )  x 

or,  x  »  (mx/mi  +  m2)  X!  (6) 

Since  both  the  conductor  and  platform  are  moving  at  the  same 
velocity#  the  value  of  x  is  the  value  of  ^  in  equation  3. 

Finally,  an  equation  is  needed  to  determine  if  the  plat¬ 
form  and  conductor  have  separated.  The  contact  force  can  be 
used  to  perform  this  determination.  The  diagram  below 
illustrates  the  contact  force, f(t). 


b2  x 
k2x2. 

m2  x 


FIGURE  1.2  CONTACT  FORCE  COMPONENTS 


7.  Contact  force: 

f(t)  b2x  +  k2x2 


m  z  x 


(7) 


-13- 


In  summary,  the  ocean  platform  has  been  modeled  and  the 
mathematical  equations  have  been  written  to  define  the  motion 
of  the  platform  system  at  any  point  in  time.  The  next 
chapter  illustrates  the  process  of  non-dimensionalizing  these 
equations  which  will  lend  them  to  analysis  on  a  more  general 
level.  A  numerical  integration  method  of  solution  will  then 
be  presented  after  which  the  results  of  the  analysis  will 


be  determined. 
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CHAPTER  2 

Non-dimensional ization  of  the  System  Equations  of  Motion 

The  removal  of  dimensions  from  the  system  equations  of 
motion  performs  two  important  functions.  First,  the  dimension¬ 
less  parameters  of  the  system  have  a  greater  application. 
Secondly,  systematic  variation  of  the  dimensionless  parameters 
will  provide  data  which  can  be  compared  to  known  standards 
when  the  excitation  is  periodic. 

Eighteen  parameters  are  contained  in  equations  1  through 
6  in  chapter  1.  The  equations  can  be  non-dimensionalized 
by  proper  selection  of  reference  quantities  for  length,  mass 
and  time.  The  author  defines  F/k  1(  mi  and  2ir/i»;01  as  the 
reference  quantities  for  length,  mass  and  time  respectively. 

The  variables  in  equations  1  through  6  for  displacement, 
velocity,  acceleration,  total  gap  distance  and  time  are  made 
dimensionless  by  division  and  differentiation  using  these 
reference  values.  The  conductor  parameters  of  mass,  viscous 
damping  coefficient  and  spring  constant  are  non-dimensional¬ 
ized  by  dividing  by  the  respective  platform  parameters.  After 
substitution  of  the  new  variables  into  equations  1,2,3  and 
7,  the  damping  cpefficient  of  the  platform  is  divided  by 
product  of  two  times  Uj  times  w01  .  This  ratio  is  defined  as 
the  damping  factor,^,  and  is  dimensionless.  The  contact 
force  variable  becomes  dimensionless  after  all  the  parameters 
of  equation  7  have  been  non-dimensionalized. 


For  clarity,  the  contact  force  variable  is  redefined  F  to 
identify  it  as  dimensionless.  The  remaining  parameter, 
excitation  frequency,  is  nondimens ional ized  by  dividing  by 
the  platform  natural  frequency.  The  dimensionless  parameters 
are  summarized  in  Table  2.1. 


Table  2.1,  Summary  of  Non-Dimensional  Parameters 


Non-dimensional 


Dimensional 


x1/CF/k1) 

x2/(F/k,) 

t/(2ir/w0i)* 

k2/k, 

mj/mj 

b2/b, 

b  i  /2ra i a) 0  j 

w/ui0 1 

1/CF/kJ  dx, 

1  / t  F/k  x )  d2xx 
1  /  ( F/k  x )  dx2 
1/CF/kx)  d2x2 
1/ (2tt/(i)o  i)  dt 
1/  (  4  it2  /  u)  i  i  )  dt2 
F/mi  g 

(do  -  di) / (F/k  i) 


(kx/mx)51 
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1.  Platform  motion  moving  separately: 

Xi  *  (1/mj)  (F  sin  tut  -  -  ktXj)  (1) 

NON- DIMENSIONAL I  ZED 

f  *  4ir2(sin  2iryiT-  r)  -  4irCj.f  (1.1) 

2.  Conductor  motion  moving  separately: 

xj  »  (l/m2)(-b2fc2  -  k2x2)  (2) 

NON-DIMENSIONALI2ED 

c  -  (1/y)  (41^5,4  -  4i:2pc)  (2.1) 

3.  Combined  platform  /  conductor  motion: 

xi  *  (l/mi+m2)(F  sin  tut  -(bi+b2)X!  -  klxl-k2x2)  (3) 
NON-DIMENSIONALIZED 


r*  (l/l+p)  (4m2  (sin  -  r  -  pc) -4-ir(;1f+-\i>Clr))  (3.1) 

4.  Gap  distances  at  any  time.t: 

gapj(t)  =  gapi(t-l)  ♦  (x2(t)  -  x2 (t-1)  - 


(xi(t)  -  xx(t-l))  (5a) 

gap2(t)  ■  do  -  di  -  gapjft)  (5b) 

NON-DIMENSIONALIZED 

gap i (t)  =  gap i (t-1)  +  (c(t)  -  c(T-l))- 

(r(T)  -  r(T-l))  •  (5a.l) 

gap 2  (t)  ■  D.  -  gapi(x)  _ _ _ C5b.l) 

5.  Conservation  of  momentum: 

x  *  (mj/m!  +  mj  Xj  (6) 


NON-DIMENSIONALIZED 
f  *  (1/(1  )  r 


(6.1) 
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Equation  6.1  appears  inconsistent.  However,  it  should  be 
interpreted  as  the  velocity  the  moment  after  a  collision  is 
equal  to  a  constant  times  the  value  of  the  velocity  just 
before  the  collision. 


Contact  Force 

f(t)  -  b2x  ♦  k2x2  -  m2x 

(7) 

Non-dimensional ized 

f  ■  Kj*  +  pc  -  (l/4ir2)yr 

(7.1) 

The  equations  of  motion  describing  the  platfoTm/conductor 
system  have  been  non-dimensional ized  and  are  in  a  form  which 
lends  then  to  analysis  through  the  use  of  an  integrative 


method  of  solution. 
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CHAPTER  3 

Method  of  Solution  and  Error  Analysis 

Since  the  interaction  between  the  platform  and  the 
conductor  are  quite  non-linear,  analytic  methods  of  solution 
must  be  bypassed  for  numerical  means.  Of  these,  the  Runga- 
JCutta  (R-K)  is  most  fitting  for  this  analysis  because  the 
method  is  self-starting  and  easily  adaptable  to  computer 
programming . 

Classically,  the  R-K  method  expresses  the  difference 
between  the  values  of  y  at  x(n+l)  and  x(n),  where  x(n+l) 
andx(n)  are  (n+1)  and  (n)  increments  of  x  away  from  a  reference 
point  xfl,  if  the  increments  are  constant  as  they  will  be  for 
this  analysis.  Mathematically  this  statement  becomes: 

y(n+l)  -  y (n)  =  Za^  (8) 

where  the  a^'s  are  constant,  the  k.'s  are  determined  by  Taylor 
series  expansion  and  'm'  is  the  order  of  the  R-K  method  being 
used.  Cursory  examination  of  equation  8  indicates  that  higher 
orders  of  the  R-K  method  yield  a  greater  number  of  terms  for 
each  value  of  y(n+l) .  It  is  to  be  shown  that  h.gher  orders  also 
yield  greater  accuracy. 

The  principle  error  in  the  Runga-Kutta  method  is  the 
truncation  error  which  is  defined: 

T  -  c(h)m+1  ♦  o(h)m+2 

■ 

where:  m  is  the  order  of  the  R-K  method 
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h  is  the  increment  (x,  -  x0  ) 
c  is  a  constant  of  dm+1 f/dxm+1 . 

The  order  of  the  hm  +  2  term  is  generally  smaller  than  the  hm+1 
term  and  will  not  be  considered  further.  The  truncation  error 
for  any  term  may  be  defined: 

T^  «  hm**  •  (number  of  steps)  , 
where  (number  of  steps)  *  constant/h. 

Therefore , 

Ti  «  hm 

An  estimation  of  the  magnitude  of  the  truncation  error  for  a 

term  can  be  determined  by  comparison  of  the  integration  values 

using  spacing  increments  of  h  and  2h. 

T  ■  y(x.)  -  y(2hU  c(2h)m  =  c(h)m2m 
m  11 

T  *  y(x.)  -  yfh^  »  c(h)m, 
m  11 

where  (2h)  and  (h)  distinguish  between  integration  increment 
values.  Subtraction  of  the  second  equation  from  the  first 
yields: 

yiCh^  ”  yi<‘2h)  *  c(h)m(2m  -  1)  or, 

Tm  -  (yiChJ-  yi(2h))/(2m  -  1). 

After  actual  calculation,  T  =  0.0.... o,  where  b  is  the  mth 

m 

decimal  place. 

Keying  on  minimizing  truncation  error,  the  author  selected 


a  fourth  ord« r  Runga-Kutta  method  for  use  in  this  paper. 

The  numerical  approximation  equations  for  this  method  are  as 
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follows : 

y(n+l)  -  y(n)  =  |  ( k x  +  2k2  +  2k3  +  k„) 

kj  -  h  f (x(n) ,  y (n) ) 

k2  ■  h  f (x(n)  ♦  *sh,  y(n)  *  Sjkj) 

k3  *  h  f (x(n)  +  hh,  y (n)  +  *ik2) 

k*  -  h  f(x(n)  +  h,  y  (n)  +  k3) 

To  perform  the  R-K  integration  for  the  determination  of 
the  platform  and  conductor  displacements,  it  is  necessary  to 
transform  the  second  order  differential  equations  of  chapter  2 
into  first  order  differential  equations.  This  elementary 
procedure  yields  the  equations  shown  below. 


Platform  Displacement  (moving  separately) 

ri  "  r 

fx  -  r2 

•  •• 
r2  ■  r 

r2  *  4ir2(  sin  2ir  t  -  rt)  -  4-iT^rj 


Conductor  Displacement  (moving  separately) 
c1  «  c 

Cl  ■  c2 

•  •• 
c2  ■  c 

c2  ■  (1/u)  (4irKtc2  -  4tt2  pc  j ) 


(1.2) 

(1.3) 


(2.2) 

(2.3) 


T - - 


Combined  Platform/Conductor  Displacement 


t2  »  (1/1+u)  (4it2  ( sin  2irYjt  -  -  pcx) 

-4-rr C ^ i r 2  «■  <Kir2)) 
rx  -  r2 

Platform/Conductor  Gap 

gap“  (t)  =  Jap !  (x-1)  +  (Cj(t)  -  c^t-1)) 

-  (Tj (T)  -  TjCT-1)) 

Contact  Force 

f  »  ♦  pcj  -  (1/ 4tt2) yr2 


(3.2) 

(3.3) 


(5a. 2) 

(7.2) 


The  computer  program  utilized  in  this  paper  was  designed 
and  written  by  the  author.  A  description,  user's  guide  and  a 
printout  of  the  program  with  output  are  presented  in  Appendix 
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CHAPTER  4 
Analysis  Results 

As  it  was  mentioned  in  chapter  1,  the  problem  has 
been  constrained  to  maintain  manageability.  As  the  first 
iteration  toward  understanding  the  response  characteristics 
of  the  platform  system,  the  author  set  the  dimensionless 
ratios  of  the  conductor  to  platform  mass,  viscous  damping 
coefficient  and  spring  constant  to  constant  values.  The 
damping  factor,  and  the  total  gap  distance  were  also 

made  equal  to  constants.  The  excitation  force  parameter,  F01, 
was  given  three  discrete  values  for  the  purpose  described  in 
chapter  1.  The  ratio  of  excitation  frequency  to  platform 
natural  frequency,  y x ,  was  given  a  range  of  seven  values  from 
O.S  to  S.O.  This  range  yields  a  summary  plot  of  the  platform 
displacement  which  is  similar  in  appearance  to  a  standard 
transfer  function  plot.  The  values  of  the  platform  mass 
and  natural  frequency  were  estimated  from  previous  platform 
analysis  experience.  The  damping  factor  was  also  determined 
in  this  manner.  With  the  values  for  natural  frequency  and 
mass,  the  platform  spring  constant  was  calculated.  The  ratios 
between  conductor  and  platform  parameters  were  then  estimated. 
Using  these  percentages,  the  values  of  conductor  mass, 
viscous  damping  coefficient  and  spring  constant  were  calculated. 
The  values  of  the  system  parameters  are  listed  in  Table  4.1. 


.  ■  . . . . . 
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TABLE  4.1,  System  Parameter  Values 


Parameter 

Value 

* 

0.009998 

C 

0.03 

V 

0.1 

r: 

0.5, 0.8, 1.0, 1.5, 2. 0,3. 0,5.0 

p 

0.01 

F01 

0.01,0.001,0.0001 

a. 

2.07  x  104  (lb-sec2) /inch 

®2 

2.07  x  103  (lb-sec2) /inch 

3.12  x  104  ( lb-sec) /ft 

b2 

3.124  x  103  (lb-sec)/ ft 

9.08  x  104  lb/in 

^2 

9.08  x  102  lb/in 

«ox 

2.09  radians/ second 

^0  2 

0.628  radians/second 

do 

19.5  inches 

di 

17.2  inches 

The  system  was  analyzed  for  the  platform  response  over  the 
range  of  for  each  excitation  force. 

The  output  of  the  computer  program  was  two  plots. 

Each  plot  was  coded  on  the  abscissa  to  denote  the  type  of 
response  that  was  plotted  and  the  specific  values  of  Yi  and 
F01  used  for  that  computer  run.  A  complete  description  of  the 
coding  and  the  individual  displacement  vs.  T  plots  for  the 
twenty-one  parameter  combinations  are  presented  with  their 
analyses  in  Appendix  B.  Each  analysis  discusses  the  effect 
of  the  excitation  frequency  on  the  response  frequency,  the 
dynamic  equilibrium  position  either  positive  or  negative  of 
the  static  equilibrium  zero  position,  the  steady-state 
magnitude  of  the  response  and  the  primary  response  frequen¬ 
cies. 

As  an  example,  if  yt  and  F01  were  set  to  values  of 
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0.5  and  0.001  respectively,  the  resulting  platform  displace¬ 
ment  and  fourier  transform  are  shown  in  Figures  4.0  and  4.01. 
The  platform  response  has  a  period  of  2x  which  is  the  period 
of  the  excitation  frequency,  1  radian/second.  The  dynamic 
equilibrium  position  is  slightly  negative.  The  Fourier 
transform  yields  two  frequencies,  the  excitation  frequency 
and  the  platform  natural  frequency. 


Figure  4.1,  Platform  Displacement  vj.  Yj 
In  Figure  4.1  the  dimensionless  platform  displacement 
is  plotted  against  the  ratio  of  excitation  frequency  to 
platform  natural  frequency.  v(.  The  non-linearity  of 
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system  is  demonstrated  by  the  reduction  in  magnitude  as  the 
excitation  force  is  decreased  from  0.01  to  0,001,  The 
decrease  has  also  produced  a  frequency  shift  in  the  platform 
response  to  the  right.  The  coincidental  displacement  plots  for 
the  excitation  force  equal  to  0.001  and  0.0001  imply  that 
the  system  is  linear  for  these  load  levels.  Thus,  there  is  a 
breakpoint  in  the  excitation  force  level  below  which  the 
non-linear  influence  of  the  conductor  on  the  platform  natural 
frequency  is  negligible.  However,  the  non-linearity  has 
reduced  the  magnitude  of  the  response. 

On  the  following  pages.  Figures  4.2  through  4.8  present 
plots  of  the  absolute  value  of  the  Fourier  transform  of  the 
platform  displacement  against  the  ratio  of  response  frequency 
to  the  platform  natural  frequency  which  has  been  designated 
Y*.  Each  figure  displays  the  displacement  magnitudes  for  a 
specific  value  of  the  ratio  of  excitation  frequency  to  platform 
natural  frequency,  y^  .  The  congruency  of  the  plots  for  F01 
equal  to  0.001  and  0.0001  over  all  values  of  y^  reinforces 
the  statements  made  previously  concerning  the  non-linear 
effect  of  the  conductor. 
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Figure  4.5,  |R(w)|  vs.  y*;  yt»  1.5 


. 


5 


Principally,  the  analysis  has  shown  that  for  large 
values  of  excitation  force  the  system  is  non-linear.  This 
is  clearly  displayed  by  the  shift  in  the  platform  natural 
frequency.  ‘f  the  excitation  load  is  below  an  as  yet  undeter¬ 
mined  level,  the  natural  frequency  of  the  platform  is  unaffect¬ 
ed  by  the  non-linear  influence  of  the  conductor.  However, 
the  non- 1 inearity  reduces  the  steady-state  platform  magnitude. 
The  analysis  consistently  attributes  the  response  frequency 
of  the  platform  to  a  function  of  excitation  frequency  and 
the  platform  natural  frequency.  Normally,  the  platform 
would  respond  at  the  excitation  frequency  in  the  steady- 
state  mode.  However,  the  analysis  was  conducted  by  numerical 
integration  until  steady-state  was  attained.  Therefore,  when 
the  Fourier  transform  of  the  platform  displacement  was  taken, 
the  transient  effects  of  the  problem  "start-up"  were  included 
and  presented  in  the  transform  results.  This  transform  is 
more  desirable  than  a  pure  steady-state  transform  since 
actual  measurements  on  the  platform  are  taken  during  random 
excitationand  the  transient  effects  will  not  have  an  oppor¬ 
tunity  to  die  out. 

With  the  dimensionless  excitation  force  equal  to  0.01, 
the  dynamic  equilibrium  of  the  system  shifted  markedly 
negative  if  the  excitation  frequency  was  near  the  platform 
natural  frequency.  At  frequencies  above  the  platform  natural 
frequency,  the  dynamic  equilibrium  always  shifted  positively, 
but  the  shift  was  nearly  indistinguishable.  The  dynamic 
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equilibrium  position -shift  was  inconsistent  when'  the  exci¬ 
tation  frequency  equalled  the  platform  natural  frequency. 

The  period  of  the  maximum  platform  displacement  in  steady- 
state  was  coincidental  with  the  excitation  period  throughout 
the  analysis.  No  further  conclusions  can  be  drawn  from  the 
data . 

The  Runga-Kutta  numerical  integration  method  has  been 
an  effective  means  of  determining  the  dynamic  characteristics 
of  the  system.  It  is  evident  that  this  analysis  has  raised 
more  questions  than  it  has  answered  because  of  the  imposed 
parameter  constraints.  Even  with  these  limitations,  if  several 
more  of  the  parameters,  i.e.,  the  platform  damping  factor 
and  the  conductor/platform  ratios  were  varied,  the  system 
response  characteristics  could  be  more  accurately  determined. 
The  system  should  then  be  analyzed  using  an  excitation  approx¬ 
imating  the  environmental  conditions  to  which  the  platform  will 
be  subjected  to  determine  the  fundamental  flexural  baseline 
response . 

The  computer  program  has  some  minor  problems  which 
could  be  eliminated  for  greater  efficiency.  The  program 
was  designed  to  analyze  the  system  with  one  set  of  parameter 
values  per  computer  run.  However,  it  would  be  more  economical 
to  have  the  capability  to  increment  those  parameters  chosen 
for  variation  over  a  set  of  values.  This  improvement  would 
keep  the  program  on  the  computer  until  all  of  the  variations 
had  been  completed.  The  Run^a-Kutta  section  of  the  program, 
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although  correct  and  complete,  does  not  allow  for  varying 
the  incrementing  variable.  If  the  difference  between  two 
adjacent,  integrated  values  of  a  variable  reaches  a  magnitude 
that  is  too  large,  the  program  may  overflow  the  computer 
memory  capacity  and  generate  output  which  contains  only  error 
messages.  The  author  endured  this  perennial  set-back 
because  of  the  input  requirements  of  the  Fourier  transform 
and  plotting  subroutines. 
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The  computer  program  was  written  to  utilize  a  4**1  order 
Runga-Kutta  (R-K)  numerical  integration  approximation  for 
determining  the  dimensionless  displacement  of  an  ocean 
platform  with  a  centrally  located  conductor.  The  values  of  the 
displacement  are  then  printed  out,  plotted  as  a  function  of  t 
and  subjected  to  a  Fourier  transform  subroutine  whose  output 
absolute  values  are  plotted  as  a  function  of  the  u>  determined 
in  the  transform.  Because  the  equations  of  the  program  are 
non-dimensionalized ,  the  values  of  u)  used  for  the  Fourier 
transform  plot  are  a  function  of  t,  not  real  time,  t.  To 
convert  u(fourier  transform)  to  an  w(real  time) ,  one  must 
divide  the  product  of  wCfourier  transform)  times  u)(pIatform 
natural  frequency)  by  2tt.  The  program  has  been  designed  to 
compute  the  platform  response  for  one  set  of  input  parameter 
values  per  program  execution.  Thus,  if  were  assigned  a 
range  of  values,  the  user  must  run  the  program  on  the  computer 
once  for  each  value  of  yl. 

The  program  is  initiated  by  setting  the  arrays  for  the 
platform  displacement  and  velocity,  the  conductor  displacement 
and  velocity,  the  gap  values  and  time  equal  to  zero.  This 
ste >  is  followed  by  individually  setting  the  values  for  the 
dimonsionl ess  parameters.  The  values  of  the  dimensional 
parameters  for  excitation  force,  platform  spring  constant, 
platform  aperture  diameter  and  the  conductor  outside  diameter 
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are  then  entered  into  the  computer  to  calculate  the  value 
of  the  dimensionless  total  gap  distance.  The  final  initial 
conditions  are  input  when  the  individual  gap  distances  arc 
set  and  the  program  begins  execution. 

The  analysis  begins  by  computing  the  first  integration 
values  of  platform  displacement  and  velocity,  and  the 
conductor  displacement  and  velocity  using  the  'masses 
separated'  R-K  integration  equations.  Gapt  is  then  calculated 
by  subtracting  the  net  change  in  platform  displacement  from  the 
net  change  in  the  conductor  displacement  and  adding  this  result 
to  the  old  value  of  gapt.  This  new  gap;  is  subtracted  from 
the  total  gap  distance  to  determine  the  new  gap2.  Gapj  and 
gap2  are  checked  for  a  zero  or  negative  value  indicating  a 
collision  between  the  platform  and  the  conductor.  If  no 
collision  has  occurred,  time  is  incremented  and  the  integration 
takes  place  again.  If  there  has  been  a  collision,  the  value 
of  the  conductor  velocity  is  changed  to  equal  a  new  value  of 
the  platform  velocity  which  is  determined  by  conservation  of 
momentum.  The  contact  force  between  the  platform  and  the 
conductor  is  checked  for  a  zero  or  negative  value.  If  the 
contact  force  is  zero  or  negative,  the  program  returns  to  the 
'masses  separated'  R-K  integration.  If  the  contact  force  is 
greater  than  zero,  the  program  executes  the  R-K  approximation 
using  the  equations  of  motion  for  the  combined  masses.  The 
program  again  checks  for  a  zero  contact  force  and  increments 
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the  R-K  approximations  for  separated  or  combined  masses  as 
the  situation  dictates  until  the  limit  of  the  increment  count 
is  reached. 

At  this  time  the  platform  displacement  is  transferred 
into  a  plotting  subroutine  for  a  pictorial  display  of  the 
result.  The  platform  displacement  array  is  then  transformed 
into  a  complex  number  array  and  transferred  to  the  Fourier 
transform  subrout ine ( FFT) .  The  output  of  the  FFT  is  transformed 
into  an  array  of  absolute  magnitudes  and  returned  to  the 
plotting  subroutine  for  display.  Upon  completion  of  the 
plot,  the  program  ends. 
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USER'S  GUIDE 

NOTE:  This  computer  program  has  been  designed  for  execution 

on  the  INTERDATA  M70/M80  computers  of  the  M.I.T.  Joint 
Mechanical/Civil  Engineering  Computer  Facility  exclusively. 
Alterations  will  be  required  to  fulfill  the  conditions  of 
plotting  and  Fourier  transform  subroutines  at  other  facilities 

The  computer  program  is  virtually  self-contained  requir¬ 
ing  only  two  cards  for  execution  per  computer  run  once 
the  parameter  value  cards  have  been  installed.  These  cards 
provide  the  nomenclature  for  the  'x'  and  'y'  axes  of  the 
plotting  subroutine.  These  cards  are  placed  into  the  card 
deck  as  shown  below  in  the  following  order:  platform  displace 
ment  vs.  t;  absolute  value  of  platform  displacement  vs.  u . 

The  format  of  these  two  data  cards  is  2A40. 

The  platform/conductor  system  parameters  are  included  in 
the  program  main  body.  Each  parameter  may  be  varied  individ¬ 
ually  by  changing  the  appropriate  card  value. 

The  program  execution  is  controlled  by  submitting  the 

card  deck  as  follows:  //  JOB 

//  FOR 

PROGRAM  DECK  (including  subroutines 
RETURN 
END 

//  XEQ2 
*EL  FFTS 
*BC  2866 

DATA  CARDS  (TWO) 

//  END 

The  following  pages  list  the  program  and  an  example 


output . 
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COMPLEX  PWR(512) ,WTAB(256) 

INTEGER* 2  XLAB(40) 

EXTERNAL  AREA 

DIMENSION  GAP1BR(S12),QAP2BR(512),A(2,512),  XSCL(4) 

CO**»N  Rl(512) ,R2(S12) ,C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GA»WA1,PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2  KO , C2  K1 , C2  K2 , C2  K3 , L . DAMCOl , DELTAU , FO 1 
SET  INITIAL  VALUES  FOR  ARRAYS  AND  PARAMETERS 
DO  1  I*  1,512 
R1 (I)  *  0.0 
R2(I)  -  0.0 
C1(I)  -  0.0 
C2 (i)  -  n.o 
GAPIBR(I)  ?  0.0 
GAP2BR(I)  «  0.0 
A(1 , 1)  -  0.0 
A(2,I)  -  0.0 
1  TAU(I)  -  0.0 
DAMCOl  *  0.03 
L  »  0 

PI  -  3.14159625 
GAM4A1  -  1.0 
DO  -  19.5 
DELTAU  »  .1 
F01  -  .01 


FO  -  79984. 
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AK1  ■  90784. 

DI  -  17.2 
RHO  -  .01 
PSI  -  .0099998 
AMU  -  .1 
J  »  0 

DBAR  -  (DO  -  DI)/(F0/AK1) 

GAP1BR(1)  =  0.0 
GAP2BR(1)  «  DBAR 
WRITE (5, 30) 

30  £0BMIX(T1O, 'TIME' ,T24, 'PLATFORM  DISP' ,T43, 'PLATFORM  VEL' ,T60, 'COND 
2UCT0R  DISP' ,T80,' CONDUCTOR  VEL' ,T99, 'GAP1  BAR’,//) 

START  THE  PROBLEM  BY  INTEGRATING  WITH  A  4TH  ORDER  RUNGA-KUTTA  APPR 
N  -  1 

4  DO-3! I l"  N,511 

L' r 

J  ■  I  ♦  1 

TAU(J) ■ *  TAU(I)  >  DELTAU 
A(2,J)  -  TAU(J) 

CALL  RKO 
CALL  RK1 
CALL  RK2 
CALL  RK3 

R1(J)  ■  AREA  (Rl( I), DELTAU, R1K0.R1K1 ,R1K2 ,R1K3) 

A(1,J)  -  R1(J) 

R2(J)  -  AREA  (R2 (I), DELTAU, R2K0,  R2K1.R2K2.R2K3) 
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C1(J)  -  AREA  (C1(I),DELTAU,C1K0,C1K1,C1K2,C1K3) 

C2(J)  -  AREA  (  C2(I),DELTAU,C2K0,C2K1,C2K2,C2K3) 

C  CHECK  TEE  GAP  DISTANCES  BETWEEN  THE  PLATFORM  AND  THE  CONDUCTOR 
GAPIBR(J)  -  GAPIBR(I)  +(C1(J)  -  C1(I))  -  (R1(J)  -  R1(I)) 

GAP2BR(J)  -  DBAR  -  GAPIBR(J) 

IF(GAP1BR(J)  .LT.  0.0)  GO  TO  600 
700  IF(GAP1BR(J)  .EQ.  0.0)  GO  TO  100 
IF(GAP2BR(J)  .LT.  0.0)  GO  TO  300 
900  IF(GAP2BR(J)  .EQ.  0.0)  GO  TO  100 
3  CONTINUE 

100  IF (  J  .EQ.  S12)  GO  TO  400 

C  WHEN  A  GAP  GOES  TO  0.0,  CONSERVE  MOMENTUM 
N  «  L  ♦  1 
WRITE (5, 40)  TAU(N) 

40  FORMAT(/,T6,' CONTACT  OCCURRED  AT  TIME, ' ,2X,1PE12.4,/) 

IF(GAP1BR(N)  .EQ.  0.0)  GAP2BR(N)  -  DBAR 
IF(GAP2BR(N)  .EQ.  0.0)  GAPIBR(N)  -  DBAR 
R2(N)  -(1. /(I.  +AMU) )  *  R2(N) 

C2(N)  -  R2CN) 

C  INTEGRATE  THE  SYSTEM  WITH  COMBINED  MASSES 

DO  S  I-  N,S11 

C  CHECK  THE  CONTACT  FORCE  BEFORE  INTEGRATION  TO  ENSURE  NO  SEPARATION 
FI  -  2.*  PI  *  GAMMA1  *  TAU(I) 

C  ACCEL  IS  A  DUMMY  VARIABLE  FOR  THE  ACCELERATION  AT  ANY  TIME  TAU 
ACCEL  *  (4.  *  PI/ Cl .  «■  AMU))  *  (PI  *  SIN(FI) -DAMC01  *  R2(I)  -  PS 
21  •  DAMCOl  *  R2(I)  -  PI  *  RHO  *  Cl  Cl) ) 
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C  FHAT  IS  THE  VARIABLE  FOR  THE  CONTACT  FORCE 

FHAT  -  F01  *  PSI  *  DAMCOl  *  R2(I)/PI  ♦  RHO  *  F01  *  C1(I)  -  AMU  * 
2F01  *  ACCEL/  (4.  *  PI**2) 

IF(  FHAT  .EQ.  0.0)  GO  TO  IS 
IF(  FHAT  .LT.  0.0)  GO  TO  IS 
L  -  I 
J  «  L  ♦  1 

TAU(J)  *  TAU(I)  «•  DELTAU 
A(2,J)  =  TAU(J) 

CALL  RK01 
CALL  RL11 
CALL  RK21 
CALL  RK31 

R1CJ)  ■  AREA  (MCI)  .DELTAU, R1K0.R1K1  ,R1K2 ,R1K3) 

A(1,J)  -  R1(J) 

R2(J)  -  AREA  (R2(I) , DELTAU, R2 KO, R2 Kl, R2 K2 ,R2K3) 

C1(J)  -  C1(I)  ♦(R1(J)  -  Rl(I)) 

C2(J)  •  R2(J) 

GAPIBR(J)  -  GAPIBR(I) 

GAP2BR(J)  -  GAP2BR(I) 

GO  TO  20 

IS  WRITE (5, 25)  TAU(I) ,R1(I) ,R2(I) ,C1(I) ,C2(I) ,GAP1BR(I) 

25  F0RMAT(T6,1PE12,4,T25,1PE12.4,T43,1PE12.4,T61,1PE12.4,T81 , 1FE12.4, 
2T97.1PE12.4) 

IF(J  .EQ.  S12)  GO  TO  400 


V,  it* 


I 
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WRITE (5, 11)  TAU(J) 

11  FORMAT C/.T6, 'THE  CONTACT  FORCE  IS  ZERO; MASS  SEPARATION  OCCURRED 
2AT  TIME, '2X,1PE12.4,/) 

N  ■  L  ♦  1 
GO  TO  4 

20  WRITE (5, 10)  TAU(J),R1(J),R2(J),C1(J),C2(J),GAP1BR(J) 

10  F0RMAT(T6,1PE12.4,T25,1PE12.4,T43,1PE12.4,T61,1PE12.4,T81,1PE12.4, 
2T97.1PE12.4) 

5  CONTINUE 
40C  M  ■  0 
GO  TO  53 

54  DO  50  I  *  1,512 

50  PWR(I)  •  CMPLX(R1(I),0.0) 

NBITS  -  9 

NEWTB  ■  1 
INV  ■  1 

CALL  FFT  (PWR, NBITS, INV.WTAB, NEWTB) 

C  THE  VARIABLE  'S’  IS  THE  INCREMENTAL  FREQUENCY  OVER  WHICH  THE 
C  FOURIER  TRANSFORM  WILL  BE  PLOTTED,  AND  THE  CONSTANT  IN  THE  'S' 

C  EQUATION  IS  EQUAL  TO  '2*PI/(512.*DELTAU)  ' 

S  -  0.0 

DO  51  I  «  1,512 
A(1,I)  •  CABSCPWR(I)) 

A(2,I)  -  S 

51  S  •  S  ♦  .6136 
S3  NROW  »  2 
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XSCL(l)  -  0.0 
XSCL(2)  »  0.0 
XSCLC3)  -  0.0 
XSCL(4)  -0.0 
NVARS  -  1 
NPTS  =  S12 
NX  -  2 
MOVE  -  0 
LABEL  *  4 
ISCL  -  1 
FTIME  -0.0 
LOOK  -  1 

READ(8,45)  (XLAB(I),  I  =  1,40) 

45  FORMAT  (40A2) 

CALL  P ICTR ( A , NROW , XLAB , XSCL , NVARS , NPTS , NX , MOVE , LABEL , ISCL, FTIME, 
2  LOOK) 

M  «  M  ♦  1 

IF(M  .EQ.  1)  GO  TO  54 
GO  TO  35 

600  GAPIBR(J)  =0.0 
GAP2BR(J)  «  DBAR 
Cl(J)  =  R1(J) 

GO  TO  700 

800  GAP2BR(J)  "0-0 
GAPIBR(J)  =  DBAR 
C1(J)  ■  R1(J)  «•  DBAR 
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C  THE  RKO  SUBROUTINE  DETERMINES  THE  K1  COEFFICIENTS  DURING 

C  COMPONENT  SEPARATION 

SUBROUTINE  RKO 

COMMON  Rl(512) ,R2(512) ,C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAMMA1,PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1 ,R2K2 ,R2K3, 
3C2K0.C2K1 ,C2K2,C2K3, L, DAMCOl , FO 1 ,  DE LT AU 
PI2  »  PI  **2 

FI  -  2.*PI  *  GAMMA1  *  TAU(L) 

R1KO  -  R2(L) 

R2K0  •  (4.  *PI2)  *  (SIN(FI)  -  R1(L))  -  4.  *  PI  *  DAMCOl  *  R2(L) 
C1KO  -  C2(L) 

C2K0  -  -(4.  *  PI  *  PSI  *  DAMCOl  *  C2(L)/AMU)  -  (PI2  *  4.  *  RHO  * 
2C1(L) /AMU) 

RETURN 


END 


A- 1 3 

THE  RK1  SUBROUTINE  DETERMINES  THE  K2  COEFFICIENTS  DURING 
COMPONENT  SEPARATION 
SUBROUTINE  RK1 

COMMON  R1(512),R2(512),C1(512),C2(512),TAU(S12),PI  AMU , GAfWAl , PS I , 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2K0.C2K1 ,C2K2 ,C2K3,L,DAMC01 ,F01 .DELTAU 
PI2  *  PI  **2 
B  »  DELTAU  /  2. 

T  =■  TAU(L)  ♦  B 
El  -  R1(L)  ♦  B  *  R1K0 
E2  -  R2 (L)  ♦  B  *  R2K0 
E3  ■  C1(L)  ♦  B  *  C1K0 
E4  *  C2(L)  +  B  *  C2K0 
FI  =  2.  *  PI  *  GAMMA1  *  T 
R1K1  -  E2 

R2K1  »  (4.  *  PI2)  *  (SIN(FI)  -  El)  -  4.  *  PI  *  DAMC01  *  E2 
C1K1  •  E4 

C2K1  *  -(4.  *  PI*PSI*DAMC01*E4/AMU)  -  (PI2  *  4 . *RH0*E3/AMU) 

RETURN 

END 


T 


A- 14 


C  THE  RK2  SUBROUTINE  DETERMINES  THE  K3  COEFFICIENTS  DURING 
C  COMPONENT  SEPARATION 
SUBROUTINE  RK2 

COMMON  Rl(512) ,R2(512) ,C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAMMA1 ,PSI, 
2RHO, R1K0 , R1K1 , R1K2 , R1K3 , C1K0 , C1K1 , Cl K2 , C1K3 , R2K0 , R2K1 , R2K2 , R2K3 , 
3C2K0.C2K1.C2K2 ,C2K3,L, DAMCOl , DELTAU , FO 1 


A-1S 

C  THE  RK3  SUBROUTINE  DETERMINES  THE  K4  COEFFICIENTS  DURING 

C  COMPONENT  SEPARATION 

SUBROUTINE  RK3 

COMMON  Rl(512) ,R2(512) ,C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAMMA1,PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1iC3,R2KO,R2K1,R2K2,R2K3, 
3C2  K0.C2K1 ,C2K2,C2K3,L, DAMCO 1 , DELTAU , F01 
B  =■  DELTAU 
PI2  *  PI  **2 
T  *  TAU(L)  +  B 
El  *  R1(L)  ♦  B  *  R1K2 

E2  *  R2(L)  +  B  *  R2K2 

E3  =  C1(L)  ♦  B  *  C1K2 

E4  *  C2(L)  ♦  B  *  C2K2 

FI  «  2.  *  PI  *  GAMA1  *  T 
R1K3  *  E2 

R2K3  -  (4.  *  PI2)  *  (SIN(FI)  -  El)  -  4.*PI*DAMC01*E2 
C1K3  *  E4 

C2K3  *  -(4 . *PI*PSI*DAMC01*E4/AMU)  -  (PI2*  4 . *RHO*E3/AMU) 

RETURN 


END 
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C  THE  RK01  SUBROUTINE  DETERMINES  THE  K1  COEFFICIENTS  FOR  THE 
C  COMBINED  MASSES 
SUBROUTINE  RK01 

COMMON  Rl(512) ,R2(512) ,C1(512) ,C2(S12) ,TAU(512) ,PI,AMU,GAMMA1 ,PSI, 
2RHO,R1KO,R1K1,R1K2  tRlK3,ClK0,ClKl ,C1K2 ,C1K3,R2K0,R2K1 ,R2K2 ,R2K3 , 
3C2  KO ,  C2  K1 ,  C2  K2 ,  C2  K3 ,  L ,  DAMCO 1 ,  DE  LTAU ,  FO 1 
FI  *  2.  *  PI  *  GAMMA1  *  TAU(L) 

PI2  =*  PI  **2 
R1KO  *  R2(L) 

R2K0  *  (4.*PI/(1.+AMU))  *  (PI*SIN(FI) -DAMC01*R2(L) -PSI*DAMC01*R2 
2(L)  -  PI  *  RHO  *  Cl  CL) ) 

C1KO  -  Cl  CL) 

RETURN 

END 


Cl  i-. 


-n* 
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C  THE  RK11  SUBROUTINE  DETERMINES  THE  K2  COEFFICIENTS  FOR  THE 
C  COMBINED  MASSES 

SUBROUTINE  RK11 

CONWON  Rl(512) ,R2(512) ,C1(S12) ,C2(512) ,TAU(512) ,PI,AMU,GA1#4A1,PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2K0,C2K1,C2K2,C2K3,L, DAMCOl, DELTAU, F01 
B  =  DELTAU  ♦  B 
T  »  TAU(L)  ♦  B 
FI  *  2.  *  PI  *  GAMMA1  *  T 
El  *  R1(L)  +  B  *  R1K0 
E2  «  R2(L)  ♦  B  *  R2K0 
E3  =»  C1(L)  ♦  B  *  C1KO 
R1K1  *  E2 

R2K1  -  (4.  *  PI/(1.*  AMU))*(PI*SIN(FI)  -DAMCOl *E2/FSI* DAMC01*E2 
2-PI  *  RHO  *  E3) 

C1K1  *  E3 

RETURN 

END 
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C  THE  RK21  SUBROUTINE  DETERMINES  THE  K3  COEFFICIENTS  FOR  THE 
C  COMBINED  MASSES 
SUBROUTINE  RK21 

CO&WON  R1(S12) ,R2(512) ,C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAMMA1 ,PSI, 
2RHO.R1KO ,R1K1 ,R1K2 ,R1K3,C1K0,C1K1 ,C1K2 ,C1K3,R2K0,R1K1 ,R2K2 ,R2K3, 
3C2JC0,C2Kl,C2K2,C2K3,L,DAMCOl,DELTAU,F01 
B  *  DELTAU  /  2. 

T  -  TAU(L)  ♦  B 
FI  -  2.  *  PI  *  GAJMA1  *  T 
El  »  R1(L)  ♦  B  *  R1K1 
E2  -  R2CL)  *B  *  R2K1 
E3  «  C1(L)  ♦  B  *  C1K1 
R1K2  -  E2 

R2K2  *(4.*PI/ (l.+AMU)) *(PI*SIN(FI) -DAMCQ1*E2-PSI*DAMC01*E2-PI*E1- 
2PI*  RHO  *  E3) 

C1K2  -  E3 
RETURN 


END 
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THE  RK.31  SUBROUTINE  DETERMINES  THE  K4  COEFFICIENTS  FOR  THE 
COMBINED  MASSES 
SUBROUTINE  RK31 

C0M4ONRl(512) ,R2(512) ,C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAM4A1,PSI, 
2RHO,RlKO,RlKltRlK2,RlK3,ClKO,ClKl,ClK2,ClK3,R2KO,R2Kl,R2K2,R2K3, 
3C2K0,C2K1,C2K2,C2K3,L, DAMCO 1 , DE  LTAU , FO 1 
B  »  DELTAU 
T  «  TAU(L)  +  B 
FI  -  2.  *  PI  *  GAMMA1  *  T 

El  «  R1(L)  ♦  B  *  R1K2 

E2  -  R2(L)  ♦  B  *  R2K2 

E3  *  Cl (L)  ♦  B  *  C1K2 

R1K3  *  E2 

R2K3  *(4.*PI/(l.+AMU))*(PI*SIN(FI)-DAMCOl*E2-PSI*DAMe01  *E2-PI* 
2RHO  *  E3) 

C1K3  *  E3 

RETURN 

END 


t 
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C  THE  AREA  FUNCTION  COMPUTES:  Y(N+1)  -  Y(N)  FOR  R1,R2,C1,C2 
FUNCTION  AREA(B1,B2,B3,B4,B5,B6) 

AREA  *  B1  ♦  B2/6.  *  (B3  ♦  2.*B4  +  2.*B5  B6) 

RETURN 


END 
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APPENDIX  B 


In  this  appendix  the  individual  platform  system  responses 
are  presented.  Each  plot  is  serialized  below  the  abscissa 
in  a  statement,  •  Run  Code  XYZ."  The  code  is  interpretted 
as  follows:  X  =  1 ;  the  plot  is  the  result  of  a  R-K  integra¬ 

tion  . 

Y  =  1  through  7;  denotes  the  value  of  Y!  used  in 
the  R-K  integration. 

Z  *  1  through  3;  denotes  the  value  of  F01  used 
during  the  R-K  integration. 

A  complete  description  of  the  codes  is  listed  in  Table  B.l. 

TABLE  B.l,  Computer  Run  Code  Summary 


CODE 

Yi 

F01 

CODE 

Yi 

F01 

111 

0.  S 

0.01 

152 

2.0 

0.001 

121 

0.8 

0.01 

162 

3.0 

0.001 

131 

1.0 

0.01 

172 

S.O 

0.001 

141 

1.5 

0.01 

113 

0.5 

0.0001 

151 

2.0 

0.01 

123 

0.8 

0.0001 

161 

3.0 

0.01 

133 

1.0 

0.0001 

171 

S.O 

0.01 

143 

1.5 

0.0001 

112 

0.5 

0.001 

153 

2.0 

0.0001 

122 

0.8 

0.001 

163 

3.0 

0.0001 

132 

1.0 

0.001 

173 

5.  O' 

0.0001 

142 

2.0 

0.001 

B-2 


Each  plot  is  followed  by  an  analysis  of  the  platform  response 
with  respect  to  steady-state  magnitude,  the  effect  of  the 
excitation  frequency,  the  position  of  the  dynamic  equilibrium, 
positive  or  negative,  relative  to  the  static  equilibrium  zero 
position  and  the  primary  response  frequencies  as  determined  by 
a  Fourier  transform  of  the  platform  response.  The  plots  which 

•9 

were  computed  for  F01  equal  to  0.0001  are  included  without 
an  explanation  of  the  plots.  These  plots  coincided  with  the 
platform  response  plots  for  F01  equal  to  0.001  for  all  values 
of  Yj.  Therefore  no  additional  information  could  be  obtained 
from  them. 
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FIGURE  B.l,  Platform  Response  for  F01  *  0.01;  0.5 

The  steady-state  platform  displacement  is  3.5  units. 

The  period  between  maximum  displacements  is  2t  which  corres¬ 
ponds  to  the  excitation  frequency  of  1  radian/second.  The 
dynamic  equilibrium  position  is  slightly  positive.  The  Fourier 
transform  of  the  platform  response  yields  primary  frequencies 
of  0.42  and  0.84  radians/second .  These  frequencies  correspond 
to  the  difference  between  the  excitation  frequency  and  the 
conductor  natural  frequency,  and  the  first  harmonic  of  this 
frequency. 


FIGURE  B.2,  Platform  Response  for  F01=  0.01;  Yi  =  0.8 
The  steady-state  magnitude  of  the  platform  response  is 
8.87  units.  The  period  of  the  response  is  1.2Sx  corresponding 
to  the  excitation  frequency.  The  dynamic  equilibrium  position 
is  3.84  units  negative.  The  primary  frequencies  of  the  re¬ 
sponse  are  the  excitation  frequency,  the  sum  of  the  excitation 
frequency  and  the  platform  natural  frequency,  and  an  unident¬ 
ified  combination  of  frequencies  equal  to  S.7  radians/ second . 
More  comprehensive  data  is  required  to  determine  the  exact 
components  of  this  frequency. 
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FIGURE  B.3,  Platform  Response  for  F01*  0.01;  1.0 

The  average  steady-state  response  magnitude  is  5.19  units. 
The  frequency  of  the  displacement  is  equal  to  the  excitation 
frequency.  As  the  overall  displacement  moves  in  the  negative 
direction,  the  individual  cycle  amplitudes  increase  over  a 
four  cycle  period  such  that  if  the  equilibrium  position  were 
normalized,  a  beat  frequency  would  be  prevalent.  If  this 
conclusion  is  true,  the  period  over  the  four  cycles  should  be 
equal  to  one-half  the  beat  frequency  period,  or  2tt/(o)1  -  o>2)  . 
Substituting  the  natural  frequencies  of  the  platform  and 
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the  conductor  for  Wj  and  co2  respectively,  the  period  is  4.28t. 
By  inspection  of  the  displacement  plot,  the  period  of  the 
four  cycles  is  approximately  4t.  Thus,  a  beat  frequency 
caused  by  the  interaction  of  the  platform  and  the  conductor 
is  a  plausible  explanation  of  the  platform  displacment  plot. 
The  dynamic  equilibrium  position  is  2.36  units  negative.  The 
Fourier  transform  of  the  response  indicates  a  primary 
frequency  equal  to  the  excitation  frequency  and  the  platform 
natural  frequency. 
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FIGURE  B .  4 ,  Platform  Response  for  F01»  0.01;  Yi"  1.5 
The  steady-state  response  magnitude  is  0.67  units.  The 
period  over  one  response  cycle  is  0.66t  corresponding  to  the 
excitation  frequency.  However,  the  period  between  two  adjacent 
peaks  of  equal  amplitude  is  2x.  The  dynamic  equilibrium 
position  is  slightly  positive.  From  the  Fourier  transform  of 
the  response,  the  principle  frequencies  are  the  excitation 
frequency,  the  platform  natural  frequency  and  a  frequency  equal 
to  the  difference  between  these  two  frequencies.  The  period 
of  the  third  frequency  is  2t  corresponding  to  the  predominant 


period  between  displacement  peaks. 


THCNIUNCTS'  HMCOKISI 


FIGURE  B,5,  Platform  Response  for  F01-  0,01;  yx-  2,0 
The  steady-state  response  magnitude  is  0,50  units.  The 
platform  oscillates  alternately  with  a  period  which  equates  to 
the  sum  of  the  platform  and  conductor  natural  frequencies  and 
a  period  which  is  twice  the  summation  frequency.  The  variation 
in  response  period  implies  that  the  platform  and  conductor  are 
alternately  in  and  out-of-phase.  The  dynamic  is  slightly 
positive.  The  Fourier  transform  indicates  primary  frequencies 
equal  to  the  excitation  frequency  and  the  platform  natural 
frequency. 
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FIGURE  B.7,  Platform  Response  for  F01  =  0.01;  y,»  5.0 
The  steady-state  magnitude  is  0.012  units.  Uniquely, 
the  system  settled  into  an  unstable  equilibrium  position 
before  reaching  its  final,  steady-state  equilibrium  position. 
The  system  cycles  at  a  frequency  of  2.65  radians/second 
which  equates  approximately  to  the  sum  of  the  platform  and 
the  conductor  natural  frequencies.  It  appears  that  at  this 
high  excitation  frequency  the  relative  velocities  of  the 
conductor  are  at  a  level  such  that  the  platform  and  conductor 
momenta  are  nearly  equal.  When  the  conductor  and  platform 


are  out-of-phase,  the  momentum  of  the  conductor  is  sufficient 
enough  to  cancel  the  momentum  of  the  platform  causing  a 
reduction  in  platform  motion.  When  the  two  components  are 
in  phase,  the  platform  displacement  nearly  doubles.  The 
dynamic  equilibrium  position  coincides  with  the  static 
equilibrium  position.  The  Fourier  transform  produces  frequen¬ 
cies  from  the  platform  response  equal  to  the  conductor  natural 
frequency,  a  frequency  which  is  twice  the  conductor  natural 
frequency  and  a  frequency  equal  to  the  sum  of  the  platform 
and  the  conductor  natural  frequencies. 
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FIGURE  B.8,  Platform  Response  for  F01*  0.001;  y,«  0.S 
The  steady-state  magnitude  is  1.35  units.  The  period 
of  the  response  is  2r  which  equates  to  the  excitation 
frequency.  The  dynamic  equilibrium  position  is  slightly 
negative.  From  the  Fourier  transform,  the  primary  response 
frequencies  are  the  excitation  frequency  and  the  platform 
natural  frequency. 
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FIGURE  B.9,  Platform  Response  for  F01=  0,001;  y1*  0.3 
The  average  steady-state  magnitude  is  2.6  units.  The 
individual  displacement  cycles  have  a  period  of  1.25t  which 
is  the  excitation  frequency.  There  is  a  slight  beat  frequency 
superimposed  which  has  a  period  of  10r  and  a  frequency  of 
0.2  radians/second.  The  dynamic  equilibrium  position  is 
slightly  positive.  The  Fourier  transform  yields  primary 
frequencies  equal  to  the  excitation  frequency  and  the  platform 
natural  frequency  from  the  platform  response. 
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FIGURE  B.10,  Platform  Response  for  F01  =  0,001;  Yj*  1.0 
The  average  steady-state  magnitude  is  5.3  units.  The 
platform  displacement  period  is  It  corresponding  to  the 
excitation  frequency  and  the  platform  natural  frequency. 

A  heat  frequency  is  more  porminent.  The  beat  period  is  4t 
which  is  a  frequency  of  0.52  radians/ second .  The  dynamic 
equilibrium  position  is  0.55  units  negative.  From  the 
Fourier  transform,  the  primary  frequency  is  the  excitation 
frequency.  Nothing  in  the  transform  indicates  the  component 
frequencies  which  make  up  the  beat  frequency. 
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FIGURE  B.ll,  Platform  Response  for  F01=  0.001;  1.3 

The  average  steady-state  magnitude  is  1.1  units.  '.'he 
period  of  a  displacement  cycle  is  0.67t  which  equates  to  the 
excitation  frequency.  The  beat  phenomena  has  a  period  of 
24t  and  a  frequency  of  0.08  radians/ second .  The  dynamic 
equilibrium  position  is  0.07  units  positive.  The  Fourier 
transform  of  the  platform  response  indicates  primary  frequen¬ 
cies  equal  to  the  excitation  frequency,  the  platform  natural 
frequency  and  a  frequency  of  1.87  radians/ second  which  requires 
further  analysis  for  component  identification. 
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FIGURE  B.12,  Platform  Response  for  F01»  0.001;  Yj*  2.0 
The  steady-state  displacement  magnitude  is  0.50  units. 

The  period  of  the  response  cycle  cannot  be  uniquely  determined 
as  the  period  varies  over  adjacent  cycles  as  the  displacement 
amplitude  alternately  increases  and  decreases.  However,  over 
four  cycles  the  average  period  is  0.5r  which  corresponds  to  the 
excitation  frequency.  The  period  of  the  amplitude  alteration 
is  4t  which  is  a  frequency  of  0.53  radians/ second .  The 
dynamic  equilibrium  position  is  negative  slightly.  The 
Fourier  transform  of  the  response  indicates  primary  frequen- 
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cies  equal  to  the  excitation  frequency  and  the  platform 
natural  frequency. 


•  si0i9aoaaos4O4B9P 

me  parameter  km  axe  m 


FIGURE  B.13,  Platform  Response  for  F01=  0.001;  yl*  3.0 
The  average  steady-state  magnitude  is  0.2  units.  The 
displacement  period  is  0.33t  which  is  the  excitation  frequency 
The  platform  response  exhibits  a  beat  phenomena  with  a  period 
of  24 1  and  a  frequency  of  0.08  radians/second.  The  dynamic 
equilibrium  position  is  slightly  positive.  From  the  Fourier 
transform  of  the  response,  the  primary  frequencies  are  the 
excitation  frequency  and  the  platform  natural  frequency,  It 
appears  that  the  excitation  frequency  is  sufficiently  removed 
from  the  platform  natural  frequency  that  excitation  force 
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is  an  insignificant  factor.  This  response  plot  duplicates 
the  response  plot  for  F01  equal  to  0.01. 
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FIGURE  B . 14 ,  Platform  Response  for  F01*  0.001;  Yj*  5.0 
This  platform  response  plot  coincides  exactly  with 
the  plot  for  F01  equal  to  0.01.  The  steady-state  magnitude 
is  0.012  units.  The  dynamic  equilibrium  position  coincides 
with  the  static  equilibrium  position.  The  primary  frequen¬ 
cies  from  the  Fourier  transform  of  the  response  are  the 
sum  of  the  platform  and  conductor  natural  frequencies,  the 
conductor  natural  frequency  and  a  frequency  which  is  twice 
the  conductor  natural  frequency. 


